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Abstract 

In  this  research  program  we  developed  novel  analytical  and  computational 
methods  to  estimate  the  support  of  radiating  sources  and  scatterers  from 
knowledge  of  the  corresponding  far  held  radiation  pattern  or  scattering  am¬ 
plitude.  The  focus  was  the  development  of  methods  to  estimate  the  so-called 
minimum  source  region  or  scattering  support  of  a  far  held.  The  support 
of  any  source  that  generates  the  given  far  held  must  contain  this  minimum 
source  region.  The  results  were  derived  in  the  framework  of  the  Helmholtz 
equation.  The  derived  methodology  consisted  of  two  components,  one  for 
the  estimation  of  convex  support  information,  and  a  subsequent  procedure 
to  estimate  the  nonconvex  support.  The  convex  support  information  was 
obtained  using  an  exterior  inverse  diffraction  framework.  This  allows  esti¬ 
mation  of  the  minimum  convex  source  region,  which  is  the  convex  hull  of 
the  minimum  source  region.  Nonconvex  support  information  is  extracted  in 
subsequent  steps  via  a  complementary  interior  inverse  diffraction  approach. 
This  gives  nonconvex  bounds  for  the  minimum  source  region.  In  many  prac¬ 
tical  problems,  the  latter  coincides  with  the  true  source  or  scatterer  support. 
The  derived  source  and  scatterer  support  estimation  methods  had  a  purely 
analytical  form  as  well  as  a  practical  computational  variant.  These  methods 
were  validated  with  analytical  and  computational  examples  during  the  course 
of  this  research.  Applications  of  the  derived  methods  to  imaging  were  also 
developed.  In  addition,  this  research  also  led  to  the  investigation  of  the  fun- 
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damental  limits  of  radar  target  separation  estimation,  as  measured  via  the 
Cramer-Rao  bound,  for  multiple  scattering  systems  in  which  multiple  scat¬ 
tering  is  not  negligible.  In  addition,  this  work  led  to  a  new  change  detection 
method  to  detect  unknown  targets  or  medium  changes  in  arbitrary  random 
and  complex  media.  The  method  in  question  is  based  on  the  optical  theorem 
and  was  termed  the  optical  theorem  detector.  The  method  was  found  to 
be  quite  effective,  performing  better  than  alternative  methods  such  as  the 
energy  detector. 
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1  Introduction 


The  problem  of  estimating  the  support  or  shape  of  a  scattering  target  from 
knowledge  of  its  far  zone  radiation  pattern  (or  far  field  scattering  amplitude) 
is  relevant  to  many  radar  problems  including  target  location,  detection,  clas¬ 
sification,  and  discrimination.  Target  support  information  is  also  useful  as 
a  prior  for  the  estimation  of  target  material  properties,  which  is  relevant  for 
radar  inverse  scattering  and  imaging  applications.  In  this  project  we  devel¬ 
oped  and  validated  novel  computational  methods  based  on  multipole  expan¬ 
sions  which  allow  one  to  estimate  the  support  of  an  unknown  scattering  target 
from  knowledge  of  the  far  field  radiation  pattern.  In  this  approach,  only  one 
probing  field  is  required,  and  the  computed  support  information  includes 
both  the  conventional  convex  support  as  well  as  the  possibly  nonconvex  sup¬ 
port  of  the  target.  Methods  based  on  field  inversion  (inverse  diffraction)  as 
well  as  on  source  inversion  (inverse  source  problem)  were  developed  and  val¬ 
idated.  Application  of  the  developed  target  support  estimation  approach  to 
inverse  source  and  scattering  problems  was  also  demonstrated.  The  corre¬ 
sponding  results  of  this  main  focus  area  were  published  in  two  journal  papers 
[1] ,  [2]  and  were  further  disseminated  in  several  conference  papers  and  invited 
talks  in  academia  and  a  national  laboratory. 

In  addition,  during  the  course  of  this  project  we  developed  fundamental 
Cramer- Rao  bound  characterizations  providing  insight  about  the  estimability 
of  target  separation  information.  This  problem  is  pertinent  to  the  estimation 
of  the  separation  of  two  targets  or  the  size  of  an  extended  target  from  far 
field  radar  data.  The  results  on  this  companion  problem  were  reported  in  a 
journal  paper  [3].  These  findings  provide  fundamental  performance  bounds 
on  target  support  estimability  and  complement  prior  research  on  the  related 
problem  of  target  position  estimation  developed  by  the  Pi’s  group  under 
another  project  [4]. 

In  addition,  we  also  developed  a  novel  wave  physics-based  approach  for 
the  detection  of  unknown  targets  or  medium  changes  from  active  sensor  wave- 
field  data.  This  detection  methodology  is  relevant  to  target  or  change  detec¬ 
tion  in  rather  arbitrary  unknown  backgrounds  including  random  and  complex 
media.  The  physical  foundation  for  this  novel  approach  is  a  new  theory  of  the 
optical  theorem  in  arbitrary  media,  which  was  developed  during  this  project 
within  the  full  vector  electromagnetic  formalism  in  [5].  This  theory  leads 
to  three  main  variants  of  the  optical  theorem  detector,  in  particular,  three 
alternative  indicators  or  statistics  for  change  detection:  one  based  on  real  av- 
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erage  power,  one  based  on  reactive  power  associated  to  energy  storage,  and 
another  based  on  the  total  or  apparent  power.  The  particular  form  of  this 
new  approach  for  the  context  of  one- dimensional  wave  propagation  systems, 
which  is  relevant  to  transmission  lines  and  wave  systems  that  can  be  approx¬ 
imately  modeled  in  a  single  dimension,  such  as,  e.g.,  canonical  radar,  sonar, 
and  lidar  models,  was  presented  in  [6].  This  development  provides  proof  that 
for  a  single-frequency  and  single  receiver  system  the  apparent  power  form 
of  the  optical  theorem  detector  is  equivalent  to  the  conventional  energy  de¬ 
tector.  On  the  other  hand,  the  two  approaches  differ  for  multi-dimensional 
data  such  as  in  multi-frequency  and  multi-sensor  systems.  The  application 
to  change  detection  is  formulated  and  validated  in  [7]  with  particular  em¬ 
phasis  on  the  detection  of  a  target  embedded  in  enclosed  environments  and 
in  random  media.  Further  results  pertinent  to  the  detection  of  targets  made 
of  active  media  (which  is  essential  for  active  cloaking-covered  targets)  were 
presented  in  a  recent  paper  [8].  Another  paper  is  under  preparation  which 
further  expands  the  applicability  for  transient  pulse  change  detection  of  ar¬ 
bitrary  nonlinear-media  and  time- varying  targets  [9]. 

This  report  provides  a  summary  of  the  main  results  obtained  in  this 
research  project,  along  with  references  indicating  where  these  results  can  be 
found  in  further  detail  in  the  published  literature.  We  provide  a  summary  of 
the  publication  record  associated  to  this  grant.  We  discuss  the  main  findings 
resulting  from  this  project.  We  also  provide  suggestions  for  future  research. 


2  Publications 

The  research  resulting  from  this  grant  led  to  7  journal  papers  [l]-[3],[5]-[8] 
and  6  conference  papers  [10]- [15],  including  a  paper  in  review  at  the  time  of 
submission  of  this  report  [8].  In  addition,  a  14th  paper  resulting  from  this 
grant  is  in  preparation  [9].  The  results  of  this  research  effort  were  further 
disseminated  through  invited  talks,  including  one  at  Lawrence  Livermore 
National  Laboratory  [16],  two  talks  at  AFOSR  meetings  including  one  in 
2012  [17],  and  a  number  of  talks  at  universities  [18]- [21], 
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3  Main  Findings 

The  focus  of  this  project  was  the  development  of  algorithms  for  the  estima¬ 
tion  of  the  shape  of  scatterers  from  active  wave  sensing  far  field  data.  The 
goal  was  to  demonstrate  successful  estimation  of  both  convex  and  nonconvex 
minimum  source  regions  which  are  contained  inside  the  scatterer  and  which 
in  most  inverse  problems  coincide  with  the  true  support  of  the  scatterer. 
This  main  goal  was  the  focus  of  our  research  associated  to  this  grant.  It  was 
accomplished  successfully,  and  the  main  results  were  published  in  a  num¬ 
ber  of  journal  and  conference  papers,  being  of  particular  relevance  the  two 
journal  papers  [1] ,  [2] .  In  those  papers  we  demonstrated  by  means  of  analyt¬ 
ical  multipole  theory  methods  and  computational  electromagnetic  methods 
a  new  methodology  to  estimate  the  minimum  convex  source  region  and  the 
(possibly  nonconvex)  support  of  a  scattering  target  from  knowledge  of  its 
far  scattering  amplitude.  The  results  include  many  numerical  examples  vali¬ 
dating  the  proposed  new  approach  to  estimate  the  localization  and  shape  of 
targets  from  the  far  field  data,  which  is  important  for  radar  imaging.  Ap¬ 
plications  to  imaging  were  demonstrated  in  the  papers.  The  problem  under 
consideration  and  the  main  components  of  the  formulation  and  associated 
computational  developments  are  summarized  next.  In  the  following,  we  bor¬ 
row  mostly  from  the  key  paper  [2], 

3.1  Problem  Statement  and  Relevance 

In  this  project  we  considered  the  problem  of  estimating  the  minimum  source 
region  of  a  given  far  field.  This  problem  was  formulated  in  the  framework 
of  the  scalar  Helmholtz  equation  which  captures  the  essence  of  most  elec¬ 
tromagnetic,  acoustic,  and  optical  wave  phenomena.  The  minimum  source 
region  has  the  following  meaning:  The  support  of  any  source  that  produces 
the  far  field  in  question  must  contain  this  region.  This  region  is  also  suffi¬ 
ciently  large  so  that  there  is  a  source  supported  in  it  that  generates  the  given 
far  field.  Thus  this  region  is  optimal:  it  is  the  smallest  support  of  any  source 
that  can  produce  the  far  field  in  question. 

Previous  works  addressing  this  problem  considered  the  computation  of  the 
convex  hull  of  the  minimum  source  region,  also  called  the  convex  scattering 
support  [22]  or  the  minimum  convex  source  region  [23].  Our  research  in  this 
grant  led  to  new  theory  and  algorithms  to  compute  the  actual  minimum 
source  region  including  nonconvex  minimum  source  regions  that  cannot  be 
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estimated  with  the  previous  methods.  We  also  derived  new  formulas  for 
the  minimum  source  region.  They  are  based  on  the  multipole  expansion 
[24],  They  complement  the  related  past  work  on  plane- wave-theory-based 
approach  to  estimate  the  minimum  convex  source  region  [23].  The  derived 
formulas  allow  computation  of  the  minimum  source  region  for  far  fields  that 
are  known  in  closed  form.  This  is  of  interest  in  antenna  synthesis.  It  is  also 
relevant  for  solving  certain  inverse  scattering  problems  analytically,  which  is 
a  topic  of  modern  interest  in  the  inverse  problems  community. 

The  derived  theory  and  computational  methods  for  support  estimation 
apply  to  three-dimensional  (3D)  and  two-dimensional  (2D)  wave  propagation 
systems.  The  latter  are  relevant  to  targets  that  vary  essentially  only  in  two 
transverse  coordinates  and  whose  length  in  the  other  dimension  is  much 
bigger  than  the  wavelength.  In  the  first  paper  [1]  we  considered  in  detail 
both  2D  and  3D  problems,  while  in  the  second  paper  [2]  we  emphasized 
mostly  the  2D  context.  In  both  papers  we  provided  2D  examples.  In  the 
following,  the  key  aspects  of  the  theory  are  summarized  only  for  the  2D 
context.  The  3D  theory  is  given  in  the  above  papers. 

The  principles  to  determine  the  minimum  source  region  of  a  far  field  can 
be  traced  back  to  the  work  of  Muller  ([26],  ths.  26,  27,  29)  who  used  the 
spherical  wave  or  multipole  expansion  to  show  that,  for  a  fixed  origin  O, 
there  is  a  minimum  spherical  source  region  of  radius  Rmin ,  Bmin(0 )  =  {r  G 
M3  :  r  =  |r|  <  Rm.in(0)},  such  that  in  order  for  a  source  to  produce  the  given 
far  held,  it  must  lie  in  the  interior  of  a  spherical  volume  centered  about  the 
origin  and  having  a  radius  that  is  at  least  as  large  as  Rmin ■  Furthermore, 
the  radiated  held  outside  Bmin  is  an  analytic  function  that  can  be  continued 
up  to  the  boundary  of  a  minimum  source  region  Vmin  located  inside  Bmin 
[27],  p.143.  Thus  the  far  held  can  be  inverse-diffracted  up  to  the  boundary 
of  Vmin  C  Bmin.  If  one  considers  another  origin  O'  one  can  find  another  min¬ 
imum  spherical  volume  Brnin(0'):  and  so  on.  Then  if  one  considers  a  set  of, 
say,  n  origins  and  n  associated  minimum  spherical  volumes,  the  intersection 
of  all  such  regions,  say  B„nn ,  establishes  a  sharper  bound  on  the  minimum 
source  region  Vrnin  C  Bmtn  which  lies  inside  such  intersection.  By  means 
of  this  procedure,  it  is  possible  to  uniquely  define  and  compute  the  convex 
scattering  support  or  minimum  convex  source  region,  say  Bminconv  D  Vmin, 
which  is  a  subset  of  the  convex  hull  of  any  source  radiating  the  given  far 
held.  This  minimum  convex  source  region  is  the  smallest  convex  support  of 
any  source  producing  that  held,  and  represents  important  source  localization 
information  that  can  be  inverted  uniquely  from  the  far  held  data  despite  the 
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nonuniqueness  of  the  full  inverse  source  problem  of  reconstructing  the  actual 
source  that  generated  the  given  far  field  [29,  30]. 

Previous  work  on  the  scattering  support  is  given  in  [22,  28,  31,  25,  32],  Of 
much  relevance  to  this  project  was  the  past  work  in  [28]  which  demonstrated 
that  the  convex  scattering  support  can  be  defined  via  a  Picard  test.  In 
addition,  Sylvester  [25]  demonstrated  the  existence  of  a  unique  union  of  well- 
supported  convex  sets  (UWSC  sets)  for  the  given  far  field.  A  source  of  that 
support  exists  that  generates  the  given  far  field.  Also,  this  unique  UWSC 
sets  is  a  subset  of  the  UWSC-support  of  any  source  that  produces  the  same 
far  held. 

The  computation  of  the  minimum  convex  source  region  was  also  addressed 
in  previous  work  ([23];  see  also  [27],  p.  143-146)  by  means  of  a  simpler 
methodology  based  on  the  plane  wave  expansion,  rendering  an  alternative 
approach  to  estimate  the  minimum  convex  source  region  of  a  given  far  held. 
Unlike  the  multipole  expansion  or  circular  Paley-Wiener  theorem  approaches, 
or  other  range  or  Picard  tests,  the  plane  wave  expansion  method  does  not 
depend  on  the  coordinate  system.  On  the  other  hand,  this  method  is  useful 
only  if  the  far  held  is  known  in  closed  form.  For  real  data  and  computations, 
one  needs  to  adopt  alternative  approaches  such  as  those  in  [22,  28]  and  other 
related  studies  [31,  32],  On  the  other  hand,  there  is  a  lot  of  interest  in  the 
analytical  approach  since  it  allows  handling  of  problems  where  the  far  held  is 
known  in  closed  form,  which  cannot  be  studied  in  rigor  with  computational 
approaches.  The  methods  developed  in  this  research  program  included  both 
purely  analytical  approaches  and  practical  computational  methods.  The  an¬ 
alytical  part  is  the  multipole  theory  counterpart  of  the  plane  wave  expansion 
method  in  [23] .  This  included  the  methods  to  estimate  the  convex  support  as 
well  as  the  nonconvex  support.  The  analytical  component  includes  the  mul¬ 
tipole  theory  counterpart  of  the  key  formula  (eq.  6)  in  [23],  for  estimation 
of  the  convex  support  as  well  as  the  respective  formula  for  the  nonconvex 
support.  The  practical  computational  implementation  for  both  convex  and 
nonconvex  support  estimation  was  also  developed  and  illustrated  with  exam¬ 
ples. 

It  is  important  to  comment  on  the  fundamental  and  practical  relevance 
of  minimum  source  regions.  There  are  many  useful  methods  for  imaging 
and  shape  reconstruction  of  scatterers  from  scattering  data.  The  research 
developed  thanks  to  this  grant  on  “minimum  source  regions”  is  an  important 
addition  to  the  state  of  the  art  in  this  important  area.  First  of  all,  the 
definition  of  the  minimum  source  region  is  rigorous.  It  renders  key  insight  on 
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the  information  about  the  support  of  a  source  or  scatterer  that  is  contained 
in  far  fields  of  sources  and  scatterers.  In  addition,  despite  the  nonuniqueness 
of  the  inverse  problem,  the  minimum  source  region  is  a  piece  of  information 
about  the  source  or  scatterer  that  can  be  deduced  uniquely  from  the  field 
data.  Moreover,  this  approach  has  the  particular  appeal  that  it  applies  in 
theory  to  the  extraction  of  information  about  the  support  of  a  scatterer 
using  a  single  probing  field.  Most  of  the  imaging  methods  used  in  qualitative 
inverse  scattering  apply  to  multiple  excitation  fields.  In  addition,  this  method 
performs  robustly,  even  though  it  defines  only  a  bound  for  the  true  scatterer 
support.  In  particular,  usually  the  true  support  is  similar  to  the  estimated 
minimum  source  region.  Another  aspect  that  makes  the  developed  target 
support  estimation  theory  and  algorithms  highly  relevant  is  the  fact  that  it 
is  essentially  based  on  a  Picard  test  so  that  it  can  be  extended  readily  to 
other  inverse  problems. 

3.2  Support  Estimation  Theory  and  Algorithms 

We  consider  electromagnetic  or  acoustic  radiation  and  scattering  that  can 
be  described  in  the  framework  of  the  scalar  Helmholtz  equation.  To  simplify 
the  presentation  we  consider  next  the  basic  2D  context,  however  the  general 
results  in  [1] ,  [2]  hold  for  both  2D  and  3D  problems.  We  consider  next  the 
2D  context  considered  in  [2]  which  is  relevant  to  transverse  electromagnetic 
fields.  In  this  2D  context 

{V2  +  k2)Ez(r)  =  -iuLiIz(r),  (1) 

where  Ez  represents  the  field,  Iz  is  source  distribution  which  in  the  scattering 
context  is  the  source  induced  in  the  scatterer  upon  the  given  probing  field, 
u  is  the  angular  oscillation  frequency,  p  is  the  free  space  permeability,  and  k 
is  the  wavenumber,  defined  by 


k  =  (2) 

where  e  is  the  free  space  permittivity.  In  Cartesian  coordinates  r  =  (x,  y ) 
while  in  cylindrical  coordinates  r  =  (p,  <fi). 

The  radiated  field  is  given  by  ([33],  p.598) 


E~ 


r  = 


tap 

"T 


dr'HQ^(k\r  —  r'|)A(r') 


(3) 


where  is  the  Hankel  function  of  the  first  kind  and  order  0.  From  the 
addition  theorem 

OO 

flf’Wr-r'l)  =  V  (4) 

m=— oo 

where  p<  =  min(p,  p')  and  p>  =  max(p,  p')  where  p  is  the  distance  from  r  to 
the  origin  of  coordinates  and  p'  is  the  distance  from  r'  to  the  origin. 

For  a  source  whose  support  V  is  contained  in  a  ball  of  radius  R  centered 
about  the  origin  O  we  have  in  view  of  (3,4) 

OO 

E,( r=(p,0))=-^  E  P>R  (5) 


where  the  multipole  moments 


dpi  pi  /  00- 


Using  the  large  argument  approximation  for  the  Hankel  function 


Hm(kP )  ~  \  -r-e 


i[kp—m(7v/2)—Tc/4] 


we  get  from  (14)  the  far  zone  relation 


Ez{p,  <t>) 


up  r~2 


ji(kp—  7r/4) 


4  V  tt kp 


where  the  far  held  radiation  pattern 


m=  E 


so  that  in  view  of  the  orthogonality  of  the  complex  exponentials 
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In  addition,  for  nonzero  far  fields,  there  exists  a  source  inside  the  mini¬ 
mum  source  region  Vmvn  C  V  radiating  the  given  field,  corresponding  to  the 
moments  am.  By  applying  the  above  results 

OO 

Ez(r=(p,(j>))  =  -(^  amH%\ kp)e~im<t>  p  >  Rmin  (11) 

m=— oo 

where  R,rmn  is  the  radius  of  the  smallest  ball  that  is  centered  about  the  same 
origin  and  contains  the  minimum  source  region  Vmin  C  V  of  the  given  far 
held. 

We  showed  in  the  key  papers  [1] ,  [2]  that  Rmin  is  defined  by 

2  ym 

kb) 

where  T  is  the  gamma  function  and  sup  means  supremum.  This  result  is  the 
mathematical  basis  of  the  method  developed  in  those  papers  to  estimate  the 
convex  support  of  the  source  or  scatterer.  We  also  developed  and  demon¬ 
strated  the  associated  computational  implementation  with  discrete  data.  The 
derived  methodology  to  estimate  the  convex  target  support  can  be  summa¬ 
rized  as  follows. 

The  result  (eq.(12))  holds  for  any  origin  O.  The  intersection  of  the  min¬ 
imum  spherical  source  regions  Bmin  =  {r  G  M2  :  p  <  Rmin}  f°r  different 
origins  defines  a  convex  region  that  contains  the  minimum  source  region.  Us¬ 
ing  sufficiently  many  such  origins  one  can  define  the  minimum  convex  source 
region  Bmin  conv  of  the  given  far  field.  The  method  consists  of  the  following 
steps:  1)  Consider  a  number  n  of  origins  Oa,a  =  1,2,  2)  For  each 

origin,  compute  via  the  test  (12)  the  associated  minimum  radius  Rmin(Oa ) 
and  associated  minimum  circular  region  Bmin(Oa )  =  {r  G  M2  :  |r  —  Oa\  < 
Rrnin(Oa)}.  3)  The  intersection  of  the  regions  Bmin(Oa),a  =  1,2,  ■■■,n 
defines  a  convex  region  Bmtn  bounding  the  minimum  convex  source  region 
B n/ 1 n j'onv •  i.e.,  Bmuiro7lv  G  Bmin •  The  estimate  Bmin  becomes  closer  to 
B min, conv  as  one  uses  more  sample  origins  and  their  minimum  circular  re¬ 
gions. 

Since  (11)  can  be  used  for  any  of  the  test  origins  Oa ,  it  follows  that  we 
can  uniquely  associate,  to  the  given  far  field,  a  valid  radiated  field  that  is 
well  defined  everywhere  outside  the  respective  convex  region  Bmin.  This  is 
exploited  in  the  estimation  of  the  possibly  nonconvex  support  of  the  target 
as  we  outline  in  the  following. 


am |  ^  0}.  (12) 


=  sup{5  G  M+  :  lim,n_s.oor(m) 
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If  the  origin  O  is  such  that  V  is  contained  outside  a  circle  of  radius  Rmt 
centered  at  the  origin,  then  according  to  (3,4)  the  held  in  the  interior  of  this 
circle  of  radius  Rmt  is  given  by 

OO 

Ez{r  =  {p,<j>))  =  -^  Y,  9mJm(kp)e-im *  p  <  Rint  (13) 


where 

poo  p2n 

9m  =  /  dp'p'  /  d^'H^(kp')eim*'h(p',4>’)-  (14) 

J  Rmt  Jo 

Moreover,  there  is  a  source  inside  Vmin  C  V  radiating  the  same  held  outside 
V,  i.e., 

OO 

E,(r  =  (pA))  =  ~  V  9mJm(kp)e-'m*  P  <  KS,,  (15) 

m=— oo 

where  R™fn  is  the  radius  of  the  largest  ball  centered  about  the  same  origin 
that  is  disjoint  to  the  minimum  source  region  Vmin  and  is  tangential  to  it  in 
a  way  that  defines  a  concave  boundary  of  Vmin. 

We  demonstrated  in  [1] ,  [2]  that  R™fn  is  defined  by 

/ ka\m 

R‘min  =  inf(a  G  R+  :  hmm_^oo  f  —  j  \gm\/T(m  +  1)  ±  0}.  (16) 

The  procedure  to  estimate  the  minimum  source  region  associated  to  the 
far  held  of  the  radiating  source  or  scatterer  under  consideration  consists 
of  two  parts:  First  an  estimate  of  the  convex  support  is  obtained.  This 
is  followed  by  refinement  of  the  target  support  estimate  by  means  of  the 
nonconvex  support  method  based  on  (16)  or  its  computational  counterpart 
in  the  case  of  real  and  discrete  data.  The  steps  of  the  support  reconstruction 
process  can  be  summarized  as  follows: 

1)  Apply  exterior  inverse  diffraction  and  compute  the  convex  region  Bmtn 
where  Vmin  C  Bmiri)Conv  C  B„nn .  For  sufficiently  many  computational  origins 
it  is  possible  that  Bmin  B  min  .('on  v  ■ 

2)  Consider  next  reference  origins  outside  Bmin,  say  Oa,  a  =  1,  2,  •  •  • ,  n,  and 
compute  the  held  (via  exterior  inverse  diffraction)  over  circles  centered  about 
those  origins  and  having  radius  as  large  as  possible  but  not  intersecting  the 
convex  region  B„nn .  Let  us  assume  such  computational  circles  correspond  to 
radii  ba,  with  centers  Oa. 
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nonconvex 

estimate 


Figure  1:  Successive  estimation  of  the  critical  ball  of  radius  R™ln  .  based  on 
field  information  available  from  previous  iterations.  This  method  allows  the 
estimation  of  the  minimum  source  region  Vmin  associated  to  the  far  field  of 
the  given  radiating  source  or  scatterer.  (Source:  [2].) 
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3)  Using  the  value  of  the  held  on  such  circles  compute  the  respective  R™fn(a). 

4)  Compute  the  union  Uint  of  all  the  regions  =  {r  G  R2  :  |r  —  Oa\  < 

RmL(a)-  Compute  the  complement  Uint  =  {r  6  K2  :  r  ^  Umt}. 

5)  Compute  the  intersection  of  Urnt  and  Brnin ,  which  dehnes  the  revised 
estimate  Vminjest  D  Vmin  for  the  minimum  source  region: 


Vr 


min,est 


ffim.in.  f'l  U 


int 


(17) 


6)  The  information  employed  in  the  preceding  steps  is  the  held  outside  Bmin. 
As  the  estimate  Knin,^  becomes  closer  to  Vmin  than  Bmin  we  repeat  steps 
2-4  for  computational  circles  located  outside  the  available  estimate  Vrnin.est- 
Computation  of  the  held  outside  the  available  estimate  relies  on  interior  in¬ 
verse  diffraction.  The  counterpart  of  step  5  associated  to  the  second  and 
following  iterations  is  Vmintest (current)  =  I4im,est(Previ°us)  C  Umt.  By  re¬ 
peating  this  process  (see  hgure  1)  it  is  possible  to  obtain  Vmin^est  =  Vmin. 

The  numerical  implementation  of  this  method  employs  discrete  data.  For 
a  reference  origin,  for  example,  we  compute  the  far  held  radiation  pattern 
and  the  corresponding  multipole  moments.  The  number  of  signihcant  multi- 
pole  moments  corresponds  approximately  to  the  familiar  truncation  formula 
am  ~  0  ,m  >  kR  where  R  is  the  radius  of  the  smallest  circle  that  circumcribes 
the  original  source  support  and  is  centered  at  the  origin.  The  reconstruction 
algorithm  is  implemented  computationally  as  a  Picard  test.  In  particular,  we 
require  that  the  functional  energy  of  the  reconstructed  held  exhibit  a  reason¬ 
able  finite  value.  A  methodology  resembling  the  familiar  L  curve  approach 
is  developed  and  illustrated  in  [2]  which  performs  this  task  effectively. 

Figure  2,  from  the  key  paper  [2],  illustrates  the  iterative  process  to  con¬ 
struct  the  estimate  of  the  convex  support  of  a  source  consisting  of  a  point 
source  plus  a  perfect  electric  conducting  cylinder,  as  shown  in  the  hgure. 
Figure  3,  from  [2],  illustrates  the  nonconvex  support  iterations  leading  to  a 
better  estimate  of  the  total  source  support.  Figure  4,  also  from  [2],  shows  a 
source  support  estimate  obtained  with  4  convex  bounding  regions  and  4  non¬ 
convex  bounding  regions.  It  was  found  that  the  performance  of  the  developed 
methodology  is  affected  by  the  regularization  procedure,  which  establishes  a 
tradeoff  between  accuracy  and  stability.  It  represents  the  practical  way  to 
handle  numerical  noise  issues  while  maintaining  sufficient  target  informa¬ 
tion.  After  exploring  different  regularization  alternatives  we  found  that  the 
method  in  [2]  was  quite  effective  and  rendered  estimates  consistent  with  con¬ 
ventional  imaging  considerations  discussed  in  [1] ,  [2]  and  in  conference  papers 
associated  to  this  grant. 
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Figure  2:  Minimum  circles  and  cummulative  estimated  convex  boundary  of 
the  minimum  source  region  for  4  origins:  (—50,0)  (top  left),  (50,0)  (top 
right),  (10,50)  (bottom  left),  and  (10,-50)  (bottom  right).  Also  shown  is 
the  correct  support  of  the  total  source,  consisting  of  a  point  source  and  a 
perfect  electric  conducting  cylinder.  (Source:  [2].) 


Figure  3:  Cummulative  estimated  nonconvex  support,  based  on  the  esti¬ 
mated  circular  boundary  corresponding  to  a  given  test  origin.  (Source:  [2].) 
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Figure  4:  Cummulative  estimated  nonconvex  support  based  on  8  test  origins 
for  the  convex  part  and  4  origins  for  the  nonconvex  part.  (Source:  [2].) 

3.3  Summary  and  Future  Directions 

In  summary,  in  this  main  focus  area  of  this  project,  we  developed  the  required 
wave  theory  and  computational  methodology  to  exactly  and  approximately 
estimate  the  minimum  source  region  of  a  given  far  held.  This  is  relevant 
to  the  imaging  and  localization  of  sources  and  scatterers.  In  the  scatterers’ 
case,  the  probing  involves  a  single  incident  held.  Our  developments  pro¬ 
vide  the  multipole  theory  counterpart  of  the  plane  wave  expansion  approach 
considered  before  [23].  The  proposed  methodology  complements  the  previ¬ 
ous  efforts  in  this  area  which  had  emphasized  only  the  minimum  nonconvex 
source  region.  The  method  derived  in  this  project  allows  the  formulation  of 
problems  in  which  the  far  held  is  expressed  as  a  multipole  series  expansion, 
which  cannot  be  handled  via  the  prior  method  in  [23].  This  holds  for  free 
space  but  can  be  further  generalized  to  nonhomogeneous  media  via  the  plane 
wave  and  multipole  expansion  generalizations  presented  in  [34], 

We  demonstrated  the  proposed  target  support  methods  with  both  an¬ 
alytical  and  computational  examples.  The  computational  implementation 
involves  regularization  to  make  the  data  essentially  hnite-dimensional.  This 
regularization  makes  the  support  estimates  more  conservative  and  more  sta- 
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blc  than  the  analytical,  infinite- dimensional  counterpart.  We  found  that  for 
certain  scatterers  and  under  homogeneous  plane  wave  excitation  the  analyt¬ 
ical  results  essentially  define  only  the  center  of  the  scatterer.  On  the  other 
hand,  under  near  field  excitation,  the  boundary  of  the  minimum  source  re¬ 
gion  associated  to  the  induced  source  in  the  same  scatterer  tends  to  be  better 
defined  by  means  of  this  methodology.  This  appears  to  be  due  to  the  greater 
information  content  about  the  target  that  is  available  in  the  far  zone  when 
the  target  is  probed  with  near  fields.  Due  to  reciprocity  the  same  results  are 
obtained  if  one  reverses  the  role  of  the  transmitter  and  the  receiver,  i.e.,  the 
probing  fields  may  be  due  to  far  zone  sources  while  the  scattering  fields  are 
sensed  in  the  near  zone,  with  identical  results.  In  addition,  in  this  research 
we  also  demonstrated  the  use  of  the  support  estimation  step  as  a  prior  for  in¬ 
verse  source  problems,  with  encouraging  results  (a  detailed  example  is  given 
in  [2]). 

Possible  future  directions  include  a  purely  computational  variant  of  the 
methods  derived  in  this  project,  which  employs  manifolds  of  realizable  fields, 
corresponding  to  bounded  energy  sources.  Another  interesting  area  is  extend¬ 
ing  the  results  of  this  grant  to  partially  correlated  fields  which  are  relevant 
in  optics.  Another  important  extension  of  the  results  obtained  in  this  grant 
is  the  generalization  to  nonhomogeneous  media  which  can  be  accomplished 
in  principle  via  the  generalized  expansion  methods  proposed  in  [34], 

3.4  Other  Related  Research 

This  research  program  also  led  to  the  characterization,  via  the  Cramer-Rao 
bound,  of  the  estimability  of  target  separation  information.  The  focus  was 
multiple  scattering  systems  in  which  multiple  scattering  is  not  negligible  so 
that  linearizing  approximations,  in  particular,  the  Born  approximation  lead¬ 
ing  to  the  A/2  rule  of  thumb  for  the  imaging  resolution  limit,  are  not  appli¬ 
cable.  The  estimation  theoretic  results  and  associated  computer  illustrations 
were  published  in  [3]. 

In  addition,  this  grant  led  to  a  new  change  detection  radar  method  based 
on  statistics  associated  to  physical  electromagnetic  energy,  as  governed  by 
the  optical  theorem.  This  can  be  considered  an  entirely  new  field,  with  rami¬ 
fications  to  all  forms  of  active  wave  sensing,  as  well  as  rather  arbitrary  media 
and  targets.  The  detailed  results  were  presented  in  [5], [7], [6],  in  upcoming 
papers  [8], [9],  and  in  many  conference  papers  and  invited  presentations. 

In  particular,  we  developed  the  fundamental  physical  principles  [5], [6], [9] 
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applied  detection  theoretic  formulation  [7] ,  computer  simulation  tools  includ¬ 
ing  multiple  scattering  [7], [8],  and  performance  receiver  operator  character¬ 
istic  (ROC)  curves  [7], [8]  of  a  new  physics-based  coherent  change  detection 
approach  that  is  based  on  a  statistic  that  has  the  meaning  of  physical  power. 
We  termed  this  new  approach  the  optical  theorem  detector.  It  is  based  on 
a  fundamental  wave  energy  conservation  result  called  the  optical  theorem 
and  reciprocity  relations  of  wave  phenomena.  This  radar  concept  is  relevant 
to  detect  unknown  targets  or  changes  in  an  unknown  background  medium 
and  is  thus  broadly  applicable  to  many  practical  military  and  national  se¬ 
curity  situations.  The  concept  in  question  is  particularly  useful  for  highly 
reverberating  environments  such  as  enclosed  environments  (caves,  tunnels, 
indoor  facilities,  urban  canyon,  ocean  waveguide,  and  so  on).  It  was  vali¬ 
dated  with  numerical  simulations  corresponding  to  canonical  waveguide- like 
environments  as  well  as  for  random  multiple  scattering  media.  The  computer 
simulation  results  suggest  that  this  method  performs  similar  to  the  matched 
filter  even  though  it  does  not  require  prior  knowledge  of  the  sought  target. 
The  obtained  computer  simulation  results  suggest  that  this  method  tends  to 
perform  better  than  the  conventional  approach  to  change  detection  called  the 
energy  detector.  Unlike  the  energy  detector,  which  is  based  on  mathematical 
energy,  the  optical  theorem  detector  is  based  on  physical  wave  energy  asso¬ 
ciated  to  scattering,  dissipation  in  the  scatterer,  as  well  as  energy  storage  in 
the  vicinity  of  the  scatterer  (reactive  energy). 

The  most  logical  future  direction  is  the  rigorous  determination  of  the  con¬ 
ditions  (background  medium,  probing  signal,  receiver  aperture,  etc.)  under 
which  the  optical  theorem  detector  approach  is  guaranteed  to  perform  better 
than  alternative  change  detection  approaches  such  as  the  energy  detector. 
Another  direction  is  the  possible  comparison  with  alternative  wave  based 
methods  such  as  the  time  reversal  detector.  On  the  other  hand,  the  time 
reversal  detector  uses  an  additional  probing  experiment  that  is  not  used  in 
the  optical  theorem  detector.  Thus  it  is  not  obvious  that  these  two  detec¬ 
tors  can  be  compared  directly.  Perhaps  a  normalization  of  the  data  can  be 
adopted  for  this  purpose,  e.g.,  via  the  addition  of  a  second  optical  theorem 
experiment,  which  may  be  optimized  from  the  prior  data. 
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